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ABSTRACT 

We present the first unambiguous detection of quasi-periodic wave trains within the broad pulse of 
a global EUV wave (so-called "EIT wave") occurring on the limb. These wave trains, running ahead 
o/the lateral CME front of 2-4 times slower, coherently travel to distances >i?©/2 along the solar 
surface, with initial velocities up to 1400 kms -1 decelerating to ~650kms _1 . The rapid expansion 
of the CME initiated at an elevated height of 110 Mm produces a strong downward and lateral 
compression, which may play an important role in driving the primary EUV wave and shaping its 
front forwardly inclined toward the solar surface. The waves have a dominant 2 min periodicity that 
matches the X-ray flare pulsations, suggesting a causal connection. The arrival of the leading EUV 
wave front at increasing distances produces an uninterrupted chain sequence of deflections and/or 
transverse (likely fast kink mode) oscillations of local structures, including a flux-rope coronal cavity 
and its embedded filament with delayed onsets consistent with the wave travel time at an elevated 
(by ~50%) velocity within it. This suggests that the EUV wave penetrates through a topological 
separatrix surface into the cavity, unexpected from CME caused magnetic reconfiguration. These 
observations, when taken together, provide compelling evidence of the fast-mode MHD wave nature 
of the primary (outer) fast component of a global EUV wave, running ahead of the secondary (inner) 
slow component of CME-caused restructuring. 

Subject headings: Sun: activity — Sun: corona — Sun: coronal mass ejections-Sun: flares — Sun: 
oscillations — Waves 



1. INTRODUCTION 

Waves or traveling disturbances on global scales in the 
magnetized solar atmosphere have been a subject of in- 
tensive study over the past half century, partly because 
of their possible roles in transporting energy and diag- 
nostic potential in inferring physical conditions on the 
Sun via coronal seismolo gy ([U chida 1970; Roberts et al.l 
119841 : IAschwandenll2004l : iNakariakov fe Verwichte1l2005f ). 

1.1. Global EUV ("EIT") Waves 

One of the most spectacular traveling distur- 
bances in the corona is expanding extreme ultravi- 
olet (EUV) enhancements across a large fraction of 
the solar disk, discovered by the SOHO EUV Imag- 
i ng Telescop e (EIT) and often cal l ed "E IT waves" 
([Moses et al.l 119971 : iThompson et al.l I1998D or global 
EUV waves, which we adopt in this paper (or "EUV 
waves" for short). Their counterparts were also 
i maged at othe r wav e lengths, including so ft X-rays 
(jKhan fc Aur ass 2002; Narukage et al. 2002) and ra- 
dio emission (jWarmuth et al.l [2004: White & Thompson! 
|2005[ ). Traveling dist urbances i n the chromo s phere 
were detected in Ra (jMoretonl 119601 : lUchidal 119681 : 
iBalasubramaniam et a l. 2010) and He I 10803 A lines 
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(jVrsnak et al.ll200l [Gilbert fc Holzerl 120041 ). 

Global coro nal seismology using large-scale EUV waves 
(jBallail 12007). compared with local seismology using 
loop oscillations, has been little explored, because 
of their poorly understood origin and nature. Pro- 
posed mechanisms include fast-mode magnet ohydrody- 
namic (MHD) waves as the coronal counterpart of 
chromospheric Moreton waves C 



cnromospneric Moreton waves (ITnompson et all 119991: 
Wills-Dayey_fc Thompson! ll999yW ang~ 20 00|; IWu et all 
2001; Warmuth et al. 2001; Ofman & Thompson 2002; 



Grechnev et al.l 1201 lal fbl) . 

produced current 
lines 



and coronal 
shells 



mass ejection 
or successive restruc- 



(CME) 

turinq of field lines (iDelanned 120001: iDelannee et al.l 
20081: iChen et al.l 12001120051 120111 : lAttrill et al.l 120071 : 
Dai et al.l l2oTo[ T Fast modes, for example, as the 
only MHD waves that can propagate perpendicular to 
magnetic fields, can readily explain the semi-circular 
shape of global EUV waves and their observed reflection 
from coronal holes and refra c tion around active regions 
(IThompson fc Mversl 120091: iGopalswamv et all [2 009: 
T TLi et all l2012bT Tc~" 



Schmidt fc Ofmanl 120101 : iLi et al.l I2012H loimedo eT"atl 
2012). Details of relevant observations and 



mod- 



els can be found in recent reviews (Warmuth 2007, 
|20lllVrsnak fc C hverl 120081: 1 Wills- Dave v fc Attrill l l2009 : 
iGallagher fc Lond 120111: IZhukovl 120111 : iChenTl^oTl : 

iPatsourakos fc Vourlidasll2012f ). 

Early observations of global EUV waves were lim- 
ited by instrumental capacities, particularly the low 
cadence of EIT (12-18 min) and small field of view 
(FOV) of TRACE. These constraints were alleviated 
by th e STEREO/EVY Imagers (EUVI; IWuelser et all 
l2004f ) from two vantage points at improved cadence 
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of 150 s, occasionally 75 s (e.g., IVeronig et al.l 2008, 


20101: iPatsourakos & Vourlidasl 120091: iMa et al.l 2009: 


Temmer et al.ll20] 


Ll:lMuhr et aHl201lHLon£ et al.ll2011b: 


Zhao et al.l 1201 It 


Li et al.N2Q12afl. The more recent So- 



mospheric Imaging Assembly (AIA; iLemen et al .1120 ill ) 
observes the corona at seven EUV wavelengths, covering 
a wide range of temperatures, at high cadence of 12 s 
and resolution of 1//5. These capabilities have allowed 
us to study the kinematics and therm al structures of 
EUV waves in unp recedented det ail (e.g . iLiu et alll2Q10b 
IChen fc Wull20TTt ISchriiver et alll201ll ) and continue to 
shape a better understanding of this phenomenon, as we 
will see in the example presented here. 



1.2. Introduction to This Study 

Observations to date have generally shown a single, dif- 
fuse f ront in primary EUV waves (e.g. Wil ls- Dave v et all 
120071 ). We report here the first unambiguous detection of 
multiple, quasi-periodic wave trains within a broad EUV 
wave pulse, as shown in the schematic of Figure [U These 
wave trains, comprising fronts forwardly inclined toward 
the solar surface, are repeatedly launched ahead of the 
lateral CME front with a 2 min periodicity and travel at 
velocities up to ~1400kms _1 (see Figure[Sj Section 13^2]) . 
We will show that the initial CME expansion produces 
a strong downward compression, previously unnoticed, 
which may play an important role in driving the EUV 
wave (see Section l3TT|) . 

Filament os cillat ion s (|Okamotq et al.l 
iHershaw et "all 120111: iTripathi et al.l 1200 
loop deflections ( Wills-Davev & Tho mpso 



dona 

and 



1999: 



2010: 



3V _ 

Patsourakos fc Vourlidasl 120091 : ILiu et all 
Aschwanden fe Schrijverl l2011f ) triggered by global 
EUV waves or Moreton waves have been observed at 
lower cadence or less favorable viewing angles (e.g., on 
the disk). AIA's high cadence allowed us to see, for the 
first time, an uninterrupted chain sequence of deflections 
and/or oscillations of various coronal structures above 
the limb, at increasing distances on the path of the 
global EUV wave. In particular, the wave travels across 
a topological separatrix surface into a coronal cavity 
hosting a filament at an elevated velocity and instigates 
its bodily oscillations (see Figure [6j Section EJ). 

These new observations, when taken together, point 
to the true wave nature of the primary EUV wave. We 
organize this paper as follows. After an observational 
overview of the event under study in Section [2j we focus 
in Section [3] on the global EUV wave containing quasi- 
periodic wave trains running ahead of the CME, followed 
by Section S] on successive structural oscillations on the 
wave path. In Section [5j we compare these w ave trains 
with th ose in a coronal funnel behind the CME (ILiu et al.l 
1201 lbh and flare pulsations. We conclude this paper with 
Section El followed by three Appendixes for our analysis 
techniques, the kinematics and temperature and height 
dependence of the EUV wave, and the early expansion 
of the CME and its role in wave generation. 

2. OVERVIEW OF OBSERVATIONS 

This event was observed by SDO/ AIA on 2010 Septem- 
ber 8-9 in the NOAA active region (AR) 11105 on the 
northwest limb. Around 22:40 UT (see Movies 1A and 



CME Front Expanding 
Up-/Downward, Laterally 




Global EUV Wave 
Leading Front 



Fig. 1. — Schematic synopsis of the 2010 September 8-9 global 
EUV wave event, showing quasi-periodic wave trains running ahead 
of (green) and behind (red) the CME front (blue). The former is 
part of the broader EUV wave pulse that travels through a flux- 
rope coronal cavity, instigating its kink-mode oscillations. The 
CME is initiated at an elevated height and produces a downward 
and lateral compression, driving the low-corona EUV wave. The 
wave trains behind the CME originate at the flare kernel, propagate 
along a narrow funnel, and terminate at the CME front. 



(o) G0ESJ(-ray Flux (C3.3 Flare) 
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23:00 00:00 01:00 02:00 03:00 04:00 
Time (since 2010/09/08 22:30:00) 

Fig. 2. — (a) GOES soft X-ray fluxes and (b) their time deriva- 
tives (smoothed with a 2 min box car) for the C3.3 flare, (c) Helio- 
centric positions of the CME front detected by SOHO/LASCO C3 
(plus signs) and C2 (asterisks) and AIA (blue crosses, reproduced 
from Figure [l2{h)). The inset is an enlarged view of the boxed re- 
gion during the flare impulsive phase when the CME acceleration 
and global EUV wave take place. Red broken lines are parabolic 
fits to the LASCO and AIA data labeled with fitted accelerations 
a and maximum velocities % a x. The green curve is the horizontal 
position of the nor thw ard global EUV wave front in the low corona, 
taken from Figure [TuT aV which travels faster than the early CME. 



Quasi-periodic Wave Trains Within A Global EUV Wave Traveling Through A Coronal Cavity 



3 




200 400 600 800 1000 1200200 400 600 800 1000 1200200 400 600 800 1000 1200200 400 600 800 1000 1200200 400 600 800 1000 1200 
Solar X (orcsec) Solar X (orcsec) Solor X (orcsec) Solor X (orcsec) Solor X (orcsec) 

Fig. 3. — AIA images at 171 (red, top) and 193 A (green, bottom), with original on the left (see Movies 1A and IB) and running difference 
on the right (see Movie 1C for 193 A). The dashed lines indicate radial cuts R0-R14 from the eruption center (a), azimuthal cuts A0-A3 
(b), and vertical cuts V0-V11 (f) denned in Appendix A. Note EUV wave trains ahead of the CME front and oscillatory brightenings and 
darkenings anti-correlated at the two wavelengths. The four brackets in (f) marks the FOV of the enlarged inset at 211 A (blue) showing 
the coronal cavity and polar coronal hole. 



TABLE 1 

Event Timeline (2010 September 8-9) 



22:40 AR core starts gradual bulging 

23:05 Onset of rapid CME expansion 

23:05 23:30 Flare impulsive phase 

23:10 23:22 Flare X-ray pulsations (2 min period) 

23:11/23:14 Global EUV wave generated at 110 Mm to 

the North/South of the eruption center 

23:11 23:23 Quasi-periodic wave trains launched ahead of CME 

/23:14 23:26 at 110 Mm to the North/South of the eruption center 

23:11 23:29 QFP wave trains launched from flare kernel behind CME 

23:15 02:00+ Coronal cavity /filament oscillations 



IB), the active region core starts gradual bulging, accom- 
panied by the activation of a prominence that appears 
dark in EUV. At about 23:05 UT, the bulging evolves into 
a rapid expansion and leads to a CME expulsion. Simul- 
taneously, the prominence erupts and unfolds itself, ex- 
hibi ting unwinding tw is ts (described for other events by, 
e.g., iLiu et alJl2QQ9bl laL 120113 iThompson eHiIl[20TTh . 

A C3.3 class flare occurs at 23:00 UT, according to 
the GOES 1-8 A soft X-ray flux, peaks at 23:33 UT, 
and then undergoes a slow decay for several hours (see 
Figure [2j). RHESSI had a nighttime data gap during 
the flare. We thus use the spike of the GOES flux time 
derivative to define the impulsive phase, 23:05-23:30 UT. 
It is during this phase that a variety of phenomena, in- 
cluding the rapid CM E acceleration (jZhang et al.ll200H 
iTemmer et al.l l2QQ8h and the generation of the global 



EUV wave and quasi-periodic wave trains, take place. 

By 23:10 UT, a dome shaped CME front, identified 
as a bright rim surro unding a dimming region (e.g., 
IThompson et al.ll2000af h has become evident in running 
difference images (see Figure [3]). A few minutes later, 
a comparably weak intensity front appears in the low 
corona ahead of the CME flanks to the north and south. 
This weak, broad pulse, identified as a global EUV wave, 
contains quasi-periodic wave trains and propagates about 
2-4 times faster than the lateral expansion of the CME. 
The northward EUV wave travels through a coronal cav- 
ity until it reaches the polar coronal hole boundary, insti- 
gating oscillatory brightenings. B ehind the CME, addi- 
tional quasi-periodic wave trains (jLiu et al.l l2Qllbh em- 
anate from the flare along a funnel of coronal loops and 
fade away while approaching the CME front. 

Major milestones of the event are listed in Table [TJ 
We employed space-time diagrams, as detailed in Ap- 
pendix A, to analyze various features, using four sets of 
narrow cuts (slits) shown in Figures [3] and [71 These cuts 
essentially form a user-defined coordinate system to sam- 
ple images from complementing perspectives. The erup- 
tion center is defined at a height of 49" (1" = 0.731 Mm, 
R Q = 952") above the brightest flare kernel at N20.6 on 
the limb. 

3. GLOBAL EUV WAVE CONTAINING QUASI-PERIODIC 
WAVE TRAINS 



4 



Liu et al. 



Figure 3] gives an overview of the global EUV wave seen 
in space-time diagrams from azimuthal cuts AO and A3. 
The general properties of the wave are examined in Ap- 
pendix B. Here we summarize the key result: (1) Tem- 
perature variations (Appendix Bl). In general, the EUV 
wave appears as a broad stripe of darkening at 171 A and 
brightening at 193 and 211 A, indicating prompt heating 
of plasma from 0.8 to 2.0 MK followed by restoring cool- 
ing, likely caused by adiabatic compression and subse- 
quent rarefaction. (2) Kinematics (Appendix B2). Mea- 
sured at 193 A on the lowest cut AO at height h = 38", 
the leading front of the EUV wave travels with general 
deceleration at an average velocity of 650 ± 30kms _1 
to the north and ~50% slower at 370 ± 20kms _1 to 
the south. We use such velocities o to characterize the 
global EUV wave, because the 193 A wave generally sig- 
nals the first arrival of the disturbance and no systematic 
height dependence of the wave velocity has been found. 
(3) Height dependence (Appendix B3). There is a clear 
trend of delayed arrival of the EUV wave at lower heights 
in the low corona (h < 220 Mm) near the eruption. This 
indicates a wave front forwardly inclined toward the coro- 
nal base. 

3.1. EUV Wave Generation and Decoupling from CME 

To probe the origin of the EUV wave, let us examine 
the enlarged space-time diagrams from all azimuthal cuts 
in Figure The base ratio (top row) best shows the 
broad global EUV wave and post-CME dimming, while 
the running ratio (middle row) highlights fast varying 
features, such as the expanding CME/loops and quasi- 
periodic wave trains (see Section 13^2]) . 

On the lowest cut AO, the global EUV wave is first 
detected some so = 150" = 110 Mm away to the north 
and south of the eruption center, but at different times. 
The northward wave appears at 23:11:20 UT, 6 min after 
the common onset of the fast CME expansion and flare 
impulsive phase, as described in Appendix C. By that 
time, the CME front to the north near the vertical di- 
rection has reached >210 kins" 1 . The southward wave is 
generated 3 min later at the outermost edge of some low- 
lying loops rapidly pushing downward (Figures [T2l m) 
and (n)). Within the distance so, no systematic emis- 
sion change can be related to either the EUV wave or 
CME expansion. 

On higher cuts A1-A3 (Figure EJ), while the EUV wave 
signal decreases, loop expansion due to the upwardly 
growing CME becomes increasingly evident. On cuts A2 
and A3, the EUV wave can be traced further back in dis- 
tance to near the eruption center and back in time as the 
successive onsets of the CME's rapid lateral expansion. 
This suggests that the EUV wave could be the continu- 
ation into the quiet Sun region of the same disturbance 
signaling the fa st expansion (cf., differenti al expansions 
in other events; iSchriiver et~aLll2008L 120111 ). 

This rapid CME expansion may play a key role in driv- 
ing the EUV wave, as described in Appendix C. It in- 
volves a sizable volume of roughly the active region of 
fAR = 110 Mm in radius and originates at its center of a 
height ho = 110 Mm above the photosphere. Below this 
height, the CME expansion has a significant downward 
component, while the lateral expansion is negligibly slow 
and spatially confined (see Figure [12] and Movie 1C). 



Base Ratio, Azim. Cut AO (h = 38") Cut A3 (h = 258") 




23:00 23:20 23:40 00:00 00:20 23:00 23:20 23:40 00:00 00:20 
Time (UT since 2010/09/08 22:50:00) Time (UT) 



Fig. 4. — Base ratio space-time diagrams at 171 (red), 193 
(green), and 211 A (blue) from selected azimuthal cuts AO and A3 
defined in Figure O^b). The distance is measured along each cut in 
the counter-clockwise direction from the eruption center shown in 
Figure [3ja). Distances shown here are normalized to that on the 
lowest cut AO to represent the azimuthal (latitudinal) position, a 
convention taken throughout this paper, while velocity measure- 
ments were based on the original distances. The dotted lines are 
linear fits to the global EUV wave front, labeled with fitted veloci- 
ties (in kms -1 ). The horizontal dashed lines delineate the bound- 
aries within each cut of the coronal cavity and polar coronal hole, 
best seen as dark areas in the 211 A inset of Figure [3] The cross 
and plus signs on the right indicate the azimuthal positions of the 
vertical cuts defined in Figure |3ff) and of the temporal intensity 
profiles shown in Figure [9] respectively. 

Above this height, as shown in Figures [5fc) and (d), 
the lateral CME expansion experiences an acceleration 
phase with a = (0.18 ± 0.01)-(0.58 ± 0.04) kms -2 and 
maximum velocities (150 ± 10)-(300 ± 30) kms -1 , fol- 
lowed by a comparable deceleration. This makes the 
CME flanks quickly fall behind the global EUV wave 
that is 2-4 times faster, c onsistent with recently reported 
CME- wave decouplings dPatsourakos & Vourlid as 2009; 
iKienreich et~all[200l [Cheng et al.ll2012[ ). 

The 110 Mm distance of the initial EUV wave from 
the eruption cent er is comparable to those of other 
EUV waves (e.g., iVeronk et all 120081 : iLiu et al. 2010) 
and Moreton waves or type II bursts ([Warmuth et al.l 
120041 ; iMuhr et al.ll2010f h suggesting a spatially extended, 
rather than a point-like source region. Whether the EUV 
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193 A, Azimuthal Cut AO (h = 58") Cut Al (h = 88") 



Cut A2 (h=168") 



Cut A3 (h = 258") 




Run. Ratio 



23:06 23:18 23:30 23:06 
rime (UT since 2010/09/08 23:00:00) 



23:18 23:30 
Time (UT) 



23:18 23:30 
Time (UT) 



23:18 23:30 
Time (UT) 



(i) Run 
200 



ratio, Cut V8 (s=546") (j) Cut V4 (s = -551") (k) W ° ve jo™ 1 IjTjO ^f« ] °%^ r ^ (|) 193 A run, diff. 23:16:32-23:16:08 



South 



2nd Wave Front in (i) 

Inclin. Angle (top oxis) 

Arrival Time (bottom) 





23:15 23:25 23:35 23:15 23:25 

Time (UT since 2010/09/08 23:08:00) Time (UT) 



23:16:00 23:16:40 23:17:20 
Wove Front Arrival Time (UT) 



700 800 900 1000 1100 1200 
Solor-X (orcsec) 



Fig. 5. Enlarged 193 A version of Figure 0] near the eruption from all azimuthal cuts in base ratio (top) and running ratio (middle). 
The red dotted lines in (c) and (d) are parabolic fits to the lateral CME front, with fitted accelerations (in kms - 2 ) and maximum velocities 
(in kms -1 ) listed in the order of the start time of the fit (marked by an open circle). Note the acceleration followed by deceleration, 
making the CME fall behind the global EUV wave (white dotted lines). The vertical dashed line marks the onset of the flare impulsive 
phase and the rapid CME expansion at 23:05 UT. Note in (e) and (f) the quasi-periodic wave trains in the low corona to the north and 
south of the eruption, in which the white data points, shifted to the left by 5 min, are the positions of the most pronounced wave pulses, 
fitted with piecewise l inea r functions (red dotted lines) labeled by resulting velocities (in kms -1 ). Bottom: (i) and (j) Enlarged running 
ratio version of Figure fTTI for vertical cuts V8 and V4, showing delayed arrival at lower heights of the wave trains, (k) Arrival time (bottom 
x-axis) at cut V8 of the second wave front in (i) versus height (y-axis). The blue dashed line (top x-axis) is the forward inclination angle 
from the local vertical of the reconstructed wave front, which is shown in red in the running difference image of (1). 
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wave represents an MHD wave or CME caused restruc- 
turing, such a considerable distance, in general, is re- 
quired for the velocity and volume of the expanding CME 
(as the disturbance driver) to grow, in order to build up 
sufficient compression to the ambient corona to produced 
detectable emission changes. In addition, the ass ociated 
6-9 min delay (cf., 1-2 min, IVrsnak et aflfl995h seems 
reasonable fo r a piston- driven, nonlinear wave t o steepen 
into a shock (jMannll 19951 : IVrsnak fc Lulicll2000D . because 
it is on t he same order of the predicted shock forma- 
tion time (jZic et alJl2QQ8l see, e.g., their Figure 6(a)) by 
applying relevant control parameters from our measure- 
ments. 

3.2. Quasi-periodic Wave Trains Ahead of CME 

The quasi-periodic wave trains ahead of the CME ap- 
pear as successive arc-shaped pulses in the low corona 
in running difference images (see Figures |3jh) and (i), 
Movie 1C). They are best seen in running ratio space- 
time diagrams (Figure middle) as fine-structure, nar- 
row stripes of alternating brightening and darkening, 
within the broader pulse of the global EUV wave seen 
in base ratio (top). The northward wave train is partic- 
ularly evident, including coherent pulses traveling some 
500" (>Rq/2) away. The first pulse in each direction 
is located at the leading front of the global wave, ap- 
pearing at so = 150" from the eruption center ahead of 
the CME flank, followed by subsequent pulses originating 
from about the same location. Each pulse has the same 
wavelength dependence as the overall global wave noted 
above (see Appendix Bl), indicating plasma heating. 

The sharpness of these narrow pulses allowed us to 
trace detailed kinematics better than the diffuse front of 
the global EUV wave. We selected the most pronounced 
pulse from cuts AO and Al in each direction, i.e., the sec- 
ond in the north and the first in the south. We then fitted 
their space-time positions with a piecewise linear func- 
tion, as shown in red in Figures [5{e) and (f), shifted by 
—5 min to allow the original data visible. On cut AO, the 
velocities within the first 1-2 minutes are 1420 ±110 and 
500 ± 50 km s -1 in the north and south, respectively, and 
then drop by ~50% to the average velocities of the diffuse 
wave front shown in panel (a). For the northward pulse, 
there is an additional velocity change upon its arrival at 
the coronal cavity, in which it travels at 980±90kms _1 , 
about 50% faster than the 620 ± 40kms _1 velocity out- 
side, followed by a drop back to this level upon its exit 
from the cavity's far side. Such non-monotonic veloc- 
ity variations resemble tho se in other EUV waves un- 
der dif ferent circumstances (jZhukov et al. 20091 : iLi et al.l 
l2012b[ ). A similar behavior occurs on cut Al, but we 
found no systematic height dependence of the wave ve- 
locity. 

The height distribution of the quasi-periodic wave 
trains can be seen in running ratio space-time diagrams 
from vertical cuts V8 (northward) and V4 (southward), 
as shown in Figures |5ji) and (j). These waves are promi- 
nent in the low corona h < 150" ~ 110 Mm and their sig- 
nals rapidly fall with height. At least five pulses (stripes), 
uniformly spaced in time by about 2 min (128 s) can be 
identified in each direction. These stripes have negative 
slopes, indicating delayed arrival at progressively lower 
heights and wave fronts forwardly inclined toward the so- 
lar surface, same as the overall global EUV wave noted 



above (see Appendix B3). To the south, the shallower 
slope is consistent with the lower propagation velocity 
and thus longer delays. The sudden change of slope from 
negative to positive at 23:35 UT is a possible indication 
of reflection of the wave at the coronal base. 

The arrival time of the second pulse from cut V8, for 
example, is plotted against height (y-axis) in Figure [5{k) 
and fitted with a parabolic function (solid line), indicat- 
ing a delay of 53 s over a height range of 98". As detailed 
in Appendix B3, we backtracked in time this delayed ar- 
rival to reconstruct the wave front seen in the running 
difference Movie 1C that is, however, difficult to identify 
in stills. The result is shown in red in panel (1). Its for- 
ward inclination angle from the local vertical is shown in 
blue in panel (k), which increases with decreasing height, 
ranging from 17° at h = 98" to 42° at the coronal base. 

4. SEQUENTIAL OSCILLATIONS ON WAVE PATH 

As mentioned earlier, upon the arrival of the global 
EUV wave at increasing distances, local coronal struc- 
tures are sequentially deflected, setting off transverse 
oscillations, which are best seen as sinusoidal stripes 
in enlarged base-ratio space-time diagrams (Figure EJ). 
In running-ratio diagrams (Figures [5j middle), they ap- 
pear as numerous shallow-sloped, narrow stripes super- 
imposed on the steep, broad pulse of the global E UV 
wave, as noted in an earlier event (jLiu et al.l I2010L see 
their Figure 5(d)). We found oscillations lasting up to 
6 hr in a variety of structures, including the northern 
polar coronal hole boundary, with velocity amplitudes of 
8-27kms _1 and periods ranging from 12 min in short 
loops to 80 min in large-scale loops. We expect to use 
such high resolution measurements as seismological tools 
in the future to i nfer physical conditions and damping 
mechanisms (e.g.. lOfman fc A schwanden 200 21 [Of man! 
l200l iWhite fc Verwichtel[20Tl iWang et al.ll2012ft . Here 
we focus on the best manifestation of such oscillations 
in the coronal cavity to shed light on the nature of the 
EUV wave. 

4.1. Oscillations of the Coronal Cavity 

We selected three dark filament threads (labeled 
Osc. 0,1, and 2) embedded in the coronal cavity as trac- 
ers of its oscillations seen in Figure [6] We obtained the 
space-time position s(t) at the center of each thread, 
using its width as the spatial uncertainty. s(t) was 
smoothed to approximate the slowly drifting equilibrium 
position s eq (t), which was then removed to obtain the net 
oscillatory displacement As(t) = s(t) — s eq (t), as shown 
in Figure [6fd). We fitted As(t) with a damped sine func- 
tion, As(t) = Aexp[— (t— £o)/ta] sin[27r(£— to)/P], shown 
as the red solid line, where A is the initial amplitude 
at t = to, P the period, and ta the e- folding damp- 
ing time. These filament threads, although located at 
different positions, have fitted parameters in a narrow 
range, indicating a coherence consistent with the bodily 
oscillations of the entire cavity. The average parameters 
are (r A ) = (119 ±14) min, (P) = (27.6 ± 0.1) min, and 
(A) = 3'/2 ± O'.'l = (2.3 ± 0.1) Mm, giving an initial ve- 
locity amplitude (v) = (2ttA/P) = (8.8 ± 0.4) kms" 1 . 

The onsets of oscillations within the cavity further 
away from the eruption are delayed in time. For ex- 
ample, Osc. 2 at As = 58" = 42 Mm from Osc. is 
delayed by At = 44 s. This translates to a velocity of 
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Fig. 6. — (a) and (b) 193 A base ratio space-time diagrams along 
azimuthal cuts AO and Al to the north of the eruption, showing 
an uninterrupted sequence of transverse oscillations set off by the 
arrival of the global EUV wave at increasing distances, (c) Zoomed 
view of (a) showing filament oscillations within the coronal cavity, 
fitted with damped sine functions in (d). Note in (c) the delayed 
onsets of oscillations at greater distances in agreement with the 
EUV wave travel time at the measured velocity. 

As/ At = 970 km s -1 , in agreement with the elevated ve- 
locity of the wave trains within the cavity found earlier 
(see Figure Efe)). 

As detailed in Appendix D, we estimated the energies 
of the cavity oscillations and global EUV wave to be 
10 27 and 10 28 ergs, respectively, significantly lower than 
the 10 31 ergs kinetic energy of the CME. The transverse 
oscillations of the cavity are likely fast kink modes. Its 
magnetic field was thus estimated via coronal seismology 
to be on the order of ~6 G with a pitch angle of ^70°, 
suggestive of a highly twisted flux rope. 



5. COMPARING TWO TYPES OF QUASI-PERIODIC WAVE 

TRAINS 

We now examine the quasi-periodic wave trains within 
the funnel behind the CME and flare pulsations, to help 
understand the wave trains in the global EUV wave ahead 
of the CME presented in Section [321 

5.1. QFP Wave Trains in Funnel Behind CME 

Figure [71 (top panels, see Movies 2 A and 2B) shows the 
quasi-periodic wave trains behind the CME as successive 
arc-shaped intensity variations on the 1-5% level, travel- 
ing along a funnel of coronal loops rooted at the flare ker- 
nel and curved toward the south. No significant signals 
can be detected in channels other than 171 A, suggesting 
that the responsible plasma is likely near this channel's 
peak response temperature of 0.8 MK. In space-time di- 
agrams (bottom panels), these waves are steep, recur- 
rent stripes behind the bright CME front. Parabolic fits 
(dotted lines) to a prominent wave pulse in different di- 
rections indicate large decelerations of 3-4 kms -2 , with 
initial velocities at so = 150" from the flare kernel of 
(1020 ± 130)-(1240± 160) kms -1 , which are comparable 
to those of the northward quasi-periodic waves ahead of 
the CME, but twice those of the southward waves (see 
Figure [5je)). 

We identify the wave trains in the funnel as quasi- 
periodic fast propagating (QFP) waves recently discov- 
ered by AIA dLiu et al.l[20Tol [201115 IShen fc Liul [20121) 
and confirmed as fast-mode magneto sonic waves in 3D 
MHD simulations (jOfman et alJl2011[ ). The new charac- 
teristics found here are the rapid deceleration and disap- 
pearance while approaching the CME front, which will 
be discussed in Section IQl 

5.2. Fourier Analysis of Wave Trains Behind/Ahead of 
CME and Flare Pulsations 

We now apply Fourier analysis to the quasi-periodic 
waves both ahead of and behind the CME front, and 
compare them with flare pulsations. From each of the 
space-time diagrams shown in Figures [8] (left), we ex- 
tracted a data cube from the dashed box in which the 
wave signals are prominent. We then obtained its Fourier 
power (k-uj diagram) in the domain of temporal fre- 
quency v and spatial wavenumber as shown in the 
middle. Each k-uo diagram displays a bright ridge that 
represents the dispersion relation of the wave averaged 
over the extracted duration. Within a selected region 
(dashed line) around the ridge, the (^, k) positions of 
bright pixels above a threshold (10% of the maximum), 
weighted by the power, were fitted with a linear function. 
These fits (passing through the origin) , as marked by the 
two cross signs on the sides, indicate average phase ve- 
locities ofv ph = 890 ±250, 670 ±140, -380 ± 140 kms" 1 , 
which agree with the corresponding velocities measured 
in the space-time diagrams (see Figures [5] and [7j) . Here 
the +/— signs denote the propagation direction and the 
uncertainties imply scatter due to temporal variations. 

By averaging Fourier power over wave number k within 
the selected regions, we obtained the frequency distri- 
butions of these waves (Figures [5J right). Prominent 
peaks were fitted with Gaussians to find the correspond- 
ing frequencies (periods) and their 1<t uncertainties. As 
a general trend, the QFP waves behind the CME front 
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Fig. 7. — QFP wave trains seen at 171 A in the funnel rooted at the flare kernel behind the CME front. Top: Running ratio images 
(see Movies 2 A for original and 2B for running difference). The dashed (dot-dashed) lines mark the curved cuts C0-C2 (azimuthal cuts 
A0-A3). Bottom: Running ratio space-time diagrams from selected azimuthal cut A3 and curved cuts CI and C2. Overlaid dotted lines 
are parabolic fits to a single wave front seen in the respective cuts, labeled with fitted deceleration a and initial velocity vq at a reference 
distance of so = 150" (plus sign). 



are preferentially distributed at high frequencies, with 
the highest peak (marked by an asterisk) located at 
v = 21 zb 1 mHz (P = 47 ± 2 s), while the waves ahead of 
the CME are concentrated at low frequencies, with the 
highest peaks at 128 ± 18 and 212 ± 24 s for the north- 
and southward waves, respectively. The second highest 
peak of the southward waves is at 159 ± 27 s, similar to 
the highest peak of the northward waves. We thus refer 
to both with a 2 min periodicity. Such periodicities 
can also be visually identified in the space-time domain 
on the left. 

In an earlier event (|Liu et al.l l2011bh . a 3 min pe- 
riod was found in both the QFP waves in a funnel 
and X-ray and UV flare pulsa tions (jOfman fc Suil [20061 : 
iNakariakov fc Melnikovl 1 20 09). We searched for such a 
common periodicity in this event. Because of no hard 
X-ray coverage and limb occultation of UV ribbons here, 
we used the time derivative of the GOES 0.5-4 A soft X- 
ray flu x as a proxy of hard X-rays, assuming thelNeupertl 
(11968ft effect at wor k (jVeronig et al.ll2QQ5b iLiu et al.l2QQfl : 
iNing fc Caol l2010f ). As shown in Figure EJd), there 
are clear pulsations superimposed on a slower variation 
shown in red, which was then removed to give a de- 
trended light curve shown in blue. We extracted the 
detrended data within the interval defined by the dashed 
box, covering most of the impulsive phase, and obtained 
its Fourier power (panel (1)). The power is primarily 
distributed within 3-23 mHz, with concentration at low 
frequencies or long periods (2-3.5 min). This trend is 



similar to that of the EUV waves ahead of the CME (pan- 
els (j) and (k)), but differs from that of the QFP waves 
behind the CME with stronger power at high frequencies 
(panel (i)). In particular, the prominent flare period of 
143 ± 16 s agrees, within uncertainties, with the 2 min 
periodicity of the waves ahead of the CME. The wavelet 
power (panel (h)) of the flare light curve shows a peak 
of such pulsations, five of which, marked by the vertical 
dotted lines, can be identified in panel (d). On the other 
hand, despite their different general trends, the flare pul- 
sations and the waves behind the CME share some com- 
mon frequencies within their broad distributions, includ- 
ing a prominent peak at 21 zb 1 mHz or 47 ± 2 s. 

We now return to the space-time domain to further 
explore the connections between the flare pulsations and 
the two groups of quasi-periodic EUV waves. In each 
space-time diagram of Figure [8] (left), we located the on- 
set times of the first identified flare pulsation at the origin 
(s = 0, the flare kernel) and of the quasi-periodic waves 
at the distance of the near edge of the box. We then 
connected the two positions with a straight dotted line, 
followed by parallel lines corresponding to the rest of the 
five pulsations at the 143 s period. Interestingly, these 
slanted lines translate to velocities of 1460, 1350, and 
490kms _1 that agree with the initial velocities of the re- 
spective wave trains noted earlier (see Figures [5] and [7j) . 
In addition, for the northward wave train ahead of the 
CME, the slanted lines closely trace the first two pulses 
and match the onsets of the rest of the five pulses at the 



Quasi-periodic Wave Trains Within A Global EUV Wave Traveling Through A Coronal Cavity 



9 




23:12 23:18 23:24 23:30 23:36 23:42 
Time (UT since 2010/09/08 23:06:00) 



23:12 23:18 23:24 23:30 23:36 23:42 
Time (UT since 2010/09/08 23:06:00) 



10 20 30 

Frequency v (mHz) 



Fig. 8. — Fourier analysis of quasi-periodic waves and flare pulsations. Top three rows: Left: Running ratio space-time diagrams 
from curved cut CI at 171 A (a) and azimuthal cut AO at 193 A to the north (b) and south (c) of the eruption, showing waves behind and 
ahead of the CME, respectively. Middle: Fourier power {k-uj diagram) on logarithmic scales of the data cube within the boxed region on 
the left. The cross signs mark the end points of a power- weighted linear fit to the bright ridge. Right: Fourier power averaged over wave 
number k within the parallelogram in the middle. The frequencies and periods of peaks identified by red vertical lines are listed on the 
right. Selected prominent peaks are marked by asterisks and their periods are indicated by the horizontal bars on the left. Bottom row: 
(d) Time derivative of the GOES 0.5-4 A flux. Subtracting the gradual, smoothed data (red solid line) gives the detrended blue curve, 
whose wavelet power and Fourier power (within the time range defined by the dashed box) are shown in (h) and (1), respectively. The 
vertical dotted lines mark five pulsations at the dominant 143 s period. 



edge of the box. For the southward waves, such a close 
match is less evident, but the duration of these flare pul- 
sations extrapolated to the box edge overlaps with that 
of the quasi-periodic waves there. Implications of these 
coincidences will be discussed in Section 1631 

6. CONCLUDING REMARKS 

6.1. Summary 

We have presented detailed analysis of a global EUV 
wave containing quasi-periodic wave trains ahead of the 
CME front, focusing on three interrelated aspects: the 
generation of the EUV wave and its decoupling from the 
CME, sequential structural oscillations (particularly of 



the coronal cavity) triggered by the wave, and a com- 
parison of these wave trains with those behind the CME 
front and with flare pulsations. We summarize the gen- 
eral properties of various phenomena in Tableland reca- 
pitulate below our findings, the first two being reported 
for the first time: 

1. Quasi-periodic wave trains^ within the broad 
global EUV wave pulse travel ahead of the lateral 
CME front coherently to >Rq/2 along the solar 

5 Such wave trains were somewhat discernible in an earlier 
STEREO event at a lowe r cadence of 75 s, but were not analyzed 
(jPatsourakos et aLll2010al , see their Figure 11). 
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TABLE 2 

General Properties of Phenomena Involved in the Event. 

Phenomena v ( km s -1 ) Period E(ergs) 

CME lateral 150 300 vertical 990 10 5i 

Flare 2 min 

Wave Trains (WTs) ahead of CME (Global EUV Wave): 10 28 

(1) Northward Init. 1420 Ave. 650 In Cav. 980 2 min 

(2) Southward 500 370 2 min 

WTs behind CME 1020 1240 4J_s 

Cavity Osc. 9 28 min 10 2T 

surface. Their initial velocity of 1420 ± 110 kms -1 
(500 ± 50 kms -1 ) to the north (south) of the erup- 
tion decelerates to their average velocity ~50% 
slower. Their dominant 2 min period agrees with 
that of the X-ray flare pulsations. These wave 
trains are spatially confined to the low corona 
h < 110 Mm, similar to t he heights of global EUV 
waves seen by STER EO (jPatsourakos et al.ll200l 
iKienreich et al.l [2QQ9h . Their arrival is delayed at 
lower heights by up to three minutes, indicating 
wave fronts being forwardly inclined toward the 
coronal base (see Figures [5] and [5J Section 13. 2p . 
while we found no systematic height dependence 
of the wave velocity (Figure [T0|) . 

2. Sequential transverse oscillations and deflec- 
tions of local structures, upon the arrival of the 
global EUV wave at increasing distances, form an 
uninterrupted chain reaction. The best mani- 
festation is the (likely fast kink mode) oscillations 
of the coronal cavity and its embedded filament, 
with a 28 min period and 9 kms -1 velocity am- 
plitude. The delayed oscillation onsets across the 
cavity agree with the travel time of the EUV wave 
at an elevated velocity of ~980kms _1 within it. 
The estimated energies of the cavity oscillations 
and global EUV wave are 10 27 and 10 28 ergs, re- 
spectively, 10 3 -10 4 times lower than the 10 31 ergs 
kinetic energy of the CME (Figure [6j Section HJ. 

3. At low altitudes h < 110 Mm, the global EUV wave 
is first detected at so = 110 Mm to the north and 
south of the eruption center, 6-9 min after the co- 
incidental onset of the rapid CME expansion and 
flare impulsive phase. At higher altitudes, the EUV 
wave can be traced further back in distance to near 
the eruption center and back in time as the succes- 
sive onsets of rapid loop expansions involved in the 
CME (Figure [51 Section ED}. 

4. The fast expansion of the CME, with a signifi- 
cant downward component, involves a sizable vol- 
ume of about the active region of 110 Mm in radius 
with its center at an elevated height of 110 Mm 
(Figure [12j Appendix C). Negligibly slow near 
the coronal base, the lateral expansion grows with 
height, reaching 150-300 kms -1 and then deceler- 
ating, falling behind the global EUV wave of 2-4 
times faster (Figure [5]). 

5. The broad global EUV wave and its narrow quasi- 
periodic pulses appear as darkening at 171 A and 
brightening at 193 and 211 A, with the maximum 
intensity variations being generally delayed from 
193 to 211 and 171 A, suggestive of heating from 



typically 0.8 to 2.0 MK followed by restoring cool- 
ing (Figures |4] and [9j Appendix Bl). 

6. There are additional QFP wave trains (fast m ag- 
netosonic; iLiu et al.ll20TTH IQfman et al.ll20TTh of 
a shorter, dominant period of 47 ± 2 s, travel- 
ing at 1000-1200 kms -1 along the coronal funnel 
rooted at the flare kernel behind the CME front, 
near which they abruptly disappear (Figure [7J Sec- 
tion EZQ). 

6.2. Proposed Physical Picture 

We propose the following physical picture, as shown 
earlier in Figure [TJ to tie together these observations. 
After a gradual rise and bulging of the active region of 
radius tar = 150" = 110 Mm, a rapid expansion, likely 
driven by a Lorentz force, initiates near its center lo- 
cated at a height tar above the photosphere. Within 
minutes, this expansion soon involves the entire ac- 
tive region of a significant volume and becomes suffi- 
ciently fast (^OOkms -1 ) , leading to a CME expulsion 
(jPatsourakos et aLll201Qbh . An EUV wave is then gener- 
ated in the low corona at the outer edge of the expanding 
CME, so = 110 Mm from the eruption center. This dis- 
tance and the 6-9 min delay allow the CME expansion to 
grow and build up sufficient compression to the ambient 
corona or allow a nonlinear wave to steepen into a shock. 

Initially centered at a significant height of 110 Mm, the 
expanding CME volume naturally produces the sphere- 
like shape of the early EUV wave front. The elevated 
lateral expansion and underlying chromosphere of ex- 
tremely high density, essentially serving as a hard floor, 
compose an angled formation. Compared with the up- 
ward expansion into the open corona, this formation 
produces a compression to the surrounding low corona 
at an enhanced efficiency. This is similar to a pis- 
ton in a semi-enclosed volume and may play an im- 
portant role in driving the low-corona portion of the 
EUV wave that is forwardly inclined to the solar sur- 
face. The downward component of the compression^ may 
cause on-disk redsh ifts in EUV waves (jHarra et al.ll201ll 
IVeronig et al.ll2011|) and magnetic tra nsients in Moreton 
waves (jKosovichev fc Zharkoval U9 99). Such compres- 
sion can explain the preferential height distribution of 
the EUV wave in the low corona, and ca n produce den- 
sity e n hancements and adi abatic heating (jSchrijver et al.l 
120111 : iDowns et al.l [2012), which are likely responsi- 
ble for the observed emission variations, perhaps to- 
gether with other heating mechanisms. This adds to 
the importance of the lateral (vers us vertical) CME ex- 
pansion in EUV wa ve generation (jPomoell et al.l 120081 : 
iTemmer et all 120091 ). 

As the CME expands and its center rises, the ambi- 
ent corona parts to allow the CME expulsion into the 
heliosphere. Its lateral expansion at increasing heights 
then slows down and halts, partly because of the re- 
sistance of the ambience by a Lorentz restoring force. 
This produces a dimming region of a limited spatial 

6 As a possibility yet to be validated, such a downward com- 
pression co uld be physically similar to recently proposed Lorentz 
transients (Fisher et al. 2012) in flares that generate photosp heric 
magnetic field changes (Wang & Liu 2010; Sun et al. 2012) and 
even sunquakes (Zharkov et al. 2011; Hudson et al. 2012). 
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extent. However, the EUV wave proceeds at the lo- 
cal fast-mode velocity further to the rest of the corona, 
leading to the spatial decoupling from its initial driver, 
the CME. The periodicity (2 min) and coherence over 
great distances (>i? /2) of the wave trains within the 
broad global EUV wave pulse clearly indicate their true 
wave nature. Their velocities of 300-1400 kms -1 are also 
within the expected range of coronal fast-m ode speeds. 
If the simu lated "coronal Moreton waves" (jChen et al.l 
I2002L l2QQ5h repeatedly launched and running ahead of 
the CME flanks are real, our result provides the first ob- 
servational evidence of such waves. As such, global EUV 
waves can no longer be considered as single-pulse distur- 
bances as previous low-re solution observations indicated 
(jWills-Davev et aLll2007h . 

As the EUV wave propagates outward in all directions, 
the ambient coronal loops are sequentially displaced at 
increasing distances, with delays in agreement with wave 
travel times. These loops generally overshoot and start 
fast kink mode oscillations, again, driven by a Lorentz 
restoring force produced by the displaced magnetic field. 
These oscillations, signaling the first response to the ar- 
rival of the traveling disturbance, are another manifesta- 
tion of its fast-mode wave nature that is the propagation 
of a perturbation at the highest velocity supported by 
the medium. These oscillating loops do not participate 
in the CME expulsion, but are only passively perturbed, 
while maintaining their topological integrity largely un- 
changed. They ar e thus not the footprint of a CME (cf., 
lAttrill et al1l2QQ7h . 

Upon arrival at the coronal cavity, the global wave 
front and subsequent quasi-periodic wave train continue 
to travel into it and exit from its far side. Their ve- 
locities within the cavity being 50% higher than out- 
side indicate higher Alfven and fast-mode speeds, ex- 
pected from the lower plasma density an d higher mag- 
netic field strength in a flux-rope cavity (Gi bson et al.l 
12010b iBerger et af]l2Qllh . In other words, the same dis- 
turbance travels into the cavity at velocities modified ac- 
cording to the local fast-mode speed. (If the wave were 
to bypass the cavity in front of or behind it in projec- 
tion, one would not expect such a velocity change.) This 
also implies that this disturbance must have penetrated 
through a topological separatrix surface that separates 
the flux-rope, as a distinct flux system, from the sur- 
rounding corona. This contradicts the prediction that 
a CME-ca used apparent wave cannot travel acros s such 
a surface (jChen et al.l 120021 : iDelannee et~aHl2QQ7[ ). The 
EUV wave does, however, stop at the polar coronal hole 
boun dary, another separatrix s urface, as observed be- 
fore (jThompson fe Mversll2009l ). This is likely because 
of the strong reflec tion due to the large gra dient of the 
wave speeds there (jSchmidt fc OfmariBoiOh . The lack 
of reflection signature at the coronal cavity may sug- 
gest a relatively small wave speed gradient and thus a 
large transmission (versus reflection) coefficient, which 
requires further investigation. 

The QFP wave trains traveling along the coronal fun- 
nel likely origina te from the flare, as previously found 
(jLiu et al.ll2011bf ) , but with some different frequency dis- 
tributions. They are spatially confined within the CME 
bubble, which may have produced a favorable condi- 
tion for their generation and/or detection. Their sudden 
termination at the CME front is likely due to damp- 



ing generated by the discontinuity there that rapidly 
dissipates these waves, and/or due to increased dis- 
persion associated with the decrease of the fast-mode 
speed with height, as predicted in 3D MHD simulations 
(jOfman et al.ll201lh . 

6.3. Discussion 

The forward inclination of the early EUV wave front 
in the low corona h as important implications (e.g., 
iGrechnev et al.l l2011bh . An observational consequence 
is a delayed arrival by, e.g., tens of seconds to min- 
utes, of its chromospheric counterpart, a Moreton wave 
as in dicated in recent preliminary results ([White et al.l 
l201lh . While the early-stage EUV wave front is primar- 
ily shaped by its driver, as the wave travels further away, 
the propagation effect will become increasingly influen- 
tial. If the EUV wave is a true MHD wave, refraction 
can deflect its ray trajectory away from regions of high 
wave speeds toward regions of lo w wave speeds (jUchidal 
119681 : lAfanasvev fc Uralovll2Qllh . Therefore, the shape 
and temporal evolutio n of the EU V wave front has a 
seismological potential (jBallai 2007) to infer the spatial 
distribution of the fast-mode and Alfven speeds and thus 
the coronal magnetic field strength. 

The common periodicities of the quasi-periodic waves 
(especially those ahead of the CME) and flare pulsa- 
tion s an d the consistent extrapolated velocities (see Sec- 
tion !5.2p point to the possibility that a cyclic disturbance, 
likely a fast- mode M HD wave, originates from the flare 
site (cf.. lChenll2006[ ) and modulate the EUV waves both 
ahead o/and behind the CME. The different physical en- 
vironments in which these waves propagate may further 
modify their relative amplitudes by, e.g., a combination 
of phase velocity variations and wave dispersion, giving 
rise to their different observed frequency distributions. 
We also note that no clear wave signatures can be iden- 
tified within ^150" (100") from the eruption center for 
the waves ahead of (behind) the CME front. We specu- 
late this being caused by the line-of-sight integration of 
other more dominant features near the eruption and/or 
disturbance amplitudes below the instrument sensitivity 
in the early stage. Another possibility for such com- 
mon periodicities is fast kink mode oscillations of loops 
in the peripheral of the active region instigated by the 
CME eruption. In this case, the period P = 2l/(Nck) 
for the fundamental mode (N = 1) gives the loop length 
I Pc k /2 = (2 min) x (lMms~ 1 )/2 w 60 Mm, which is 
reasonable for large active region loops and comparable 
to the height range (~110 Mm, Figure |5fi)) of the low- 
corona wave trains. However, despite our best efforts, we 
could not identify such loop oscillations in AIA data. Al- 
ternatively, if these oscillations are not free, linear kink 
modes, but are driven or nonlinear, their periodicities 
would be determined by the driver or nonlinearity, not 
by the loop length. These aspects require further inves- 
tigation with 3D MHD models. 

Recent observ ati ons (jt 
IZhukov fc Aucherd 1200 ' 



Patsourakos & Vourlidas 
i Cheng et all 1201! 



arra fc Ste rling 



Kienreich et al . 
IChen fc Wul 



2009; 



2003; 



2011 



2002 



1200a 

and simulations (IChen et al .. 
120051 : iCohen et al.l l200l iDowns et al.l 1201 lL l2Q12h have 
led to a converged view about the multifaceted nature of 
global EUV waves. That is, a primary (outer) fast-mode 
MHD wave/shock, sometimes associated with an Ha 
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Moreton wave ([Thompson et al.l l2000bt lAsai et all 
l2012f h travels generally ahead of a secondary (inner) 
slow disturbance related to the CME-caused coronal 
restructuring . Along the lines of recent su ggestions (e.g., 
iDowns et al I 120111 : ISchrijver et al.l l2Qllf ) , we propose 
that this wave/non-wave hybrid mechanism is a general 
rule, rather than an exception. Under specific physical 
and observational circumstances, this general underlying 
mechanism may have different manifestations, depend- 
ing on which component has the dominant obs ervational 
signat ure. Taking the EUV wave reported by iLiu et all 
(|201Q[ ) as an example, the low-cadence, low sensitivity 
SOHO/EIT could miss the leading faint, diffuse wave 
component of 240 km s -1 and only capture the trailing 
strong, sharp non-wave fronts as slow as 80kms _1 . This 
could cause a bias in the broad vel ocity distribution of 
50-70 kms" 1 obtained from EIT ([Thompson fc Mversl 
|2009[ ). the lower end of which is incompatible with 
fast-mode waves. Likewise, if an expanding CME drives 
a fast-mode wave or shock, depending on whether the 
stand-off distance between the wave and CME front can 
be spatially resolved by the instrument, one could detect 
both or only one component. When a separatrix surface 
lies in the way, the CME-caused slow component could 
stop, while the fast-mode component would proceed 
forward. As such, the single term "EIT waves" may 
have historic ally included a variety of physically distinct 
phenomena (Wa rmuth fc Mamul201l[ ). 

To date, SDO/AIA has observed a variety of EUV 
waves, inc luding those generated by mini-CMEs, emerg- 
ing flux (Zheng et al ] I2011L 120121 ). and "coronal cy- 
clones" ( Zhang fc Liul20TTh . as well as Kelvin-Helmholtz 
instability waves produced by veloci t y she ar at CME 
flanks fsee IQfman fc Thompson! 120101, (SoTTI for the first 
AIA EUV wave ev ent reported bv ILiu et all 120101 : cf., 
iFoullon et al.l [20lTI for a later event). These new ob- 
servations, together with those presented here, add to 
the multitude of the EUV wave phenomenon in specific 
and CME/nare eruptions in general. A latest study of 
a large sample of >100 global EUV waves observed by 
both SDO/AIA and STEREO/EVVl has revealed an av- 
erage wave velocity of 700 kms -1 (N. Nitta et al. 2012, 
in preparation). This is a strong indication that the pri- 
mary EUV disturbance is a fast-mode MHD wave, a con- 
clusio n supported by many recent studies (e.g. lAsai et al.l 
|2012[ ) and particularly the quasi-periodic wave trains 
ahead of the CME discovered here. The origin of their 
periodicity will be a subject of future investigation. 

This work is supported by NASA SDO/AIA contract 
NNG04EA00C to LMSAL, WL and LO by NASA grant 
NNX11A068G, and LO by NASA grant NNX09AG10G. 
We thank Cooper Downs, Spiros Pat sour akos, P. F. 
Chen, and Keiji Hayashi for useful discussions and the 
referee for constructive comments that helped improve 
this paper. The SOHO/LASCO CME catalog is gen- 
erated and maintained at the CDAW Data Center by 
NASA and The Catholic University of America in coop- 
eration with the Naval Research Laboratory. 

APPENDIX A 
TECHNICAL DETAILS OF DATA ANALYSIS 

As shown in Figures [3] and [7J we selected four comple- 
menting sets of narrow cuts (slits) for space-time anal- 



ysis: (1) 15 radial cuts Hi (i = 0, ...,14) of 10" wide 
starting from the eruption center (h = 49" above the 
limb at N20.6 ), which serve to track the early expan- 
sion of the CME and the generation of the EUV wave; 
(2) four azimuthal cuts Ai that are 20 // -30 // wide with 
centers located at heights of h = 38", 88", 168", 258" 
and serve to track propagation at constant heights above 
the limb; (3) 12 vertical cuts Vi of 20" wide starting at 
the limb, which serve to provide height dependent arrival 
times of propagating disturbances, with cut V6 located 
at the eruption center, cut V9 within the coronal cavity, 
and cut VI 1 at the coronal hole boundary; and (4) three 
curved cuts Ci of 20" wide along the funnel in which the 
QFP waves propagate. 

Each space-time diagram was obtained by composing 
one-dimensional image profiles at different times along a 
cut. To achieve desired color scaling suitable for different 
features, we further applied base or running difference by 
subtracting an initial or immediately previous image pro- 
file, or likewise, applied base or running ratio, equivalent 
to differencing logarithmic profiles. We used mouse-clicks 
to track a moving/propagating feature whose positional 
uncertainties, depending on its sharpness, are typically 
no more than a few to ~10". 

Bear in mind that velocities measured this way are pro- 
jected velocities on the plane of the sky captured within 
the cut. They can deviate from the true plane-of-sky ve- 
locities by a factor of 1/ cos#, if the cut is at an angle 
6 from the direction of propagation normal to an inten- 
sity front. Significant deviations, however, only occur 
for large angle mis- alignments, say, a factor of 2 overesti- 
mate for = 60°. Because this event occurs on the limb, 
we expect that such projection effects are rather small. 
The selection of complementing cuts above can also help 
minimize this potential bias. 

STEREO A, 82° ahead of the Earth, provides comple- 
mentary information at a 2.5 min cadence of the CME- 
EUV wave decoupling, consistent w ith our and earlier 
results (e.g.. iPatsourakos et al.ll2009[ ). We thus chose to 
focus on new AIA observations here, but future studies 
using higher cadence STEREO data may help constrain 
the 3D geometry of quasi-periodic wave trains. 

APPENDIX B 
GENERAL PROPERTIES OF GLOBAL EUV WAVE 



Bl. 



TEMPERATURE DEPENDENCE: HEATING AND 
COOLING 



The general trend of darkening at 171 A and bright- 
ening at 193 and 211 A at the EUV wave front, followed 
by a recovery in the opposite direction, as shown 
in Figure 2J indicates heating and subsequent cooling. 
Qualitatively similar be haviors were observed by TRA CE 
and STE R EO/E UVl (|Wills-Davev fc Thompson! ^ 



Dai et all l2010h and SDO/AIA (ILiu et al.l 



1999; 



2010; 



Lon g et al.lfeoilal : iSchriiver et aTll201lD . 

We further take horizontal slices at selected positions 
from Figure [4] (plus signs there) and show the resulting 
intensity profiles in Figure [9] The large variations near 
the eruption (panels (k) and (1)) are due to the erupt- 
ing material and subsequent dimming, while the rest, 
moderate variations result from the global EUV wave. 
In particular, the CME front in panel (k) appears as a 
sharp enhancement in all three channels, suggesting an 
expected density increase. Immediately ahead of it is a 
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Fig. 9. — Base ratio temporal profiles of emission intensity from azimuthal cuts at selected positions shown in Figure 0] (plus signs there). 
The general trend of darkening at 171 A and brightening at 193 and 211 A indicates heating in the EUV wave pulse ahead of the CME. 
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Fig. 10. — (a)-(d) Positions of the north- and southward 193 A EUV wave fronts obtained from azimuthal space-time diagrams as in 
Figure |4] The sign of the southward distance is reversed. Overlaid dashed lines are parabolic (blue) and linear (red) fits, labeled with fitted 
parameters. The horizontal dotted line denotes the position where the velocity vo from parabolic fits is evaluated, (e) Height distribution 
of velocities from linear fits to wave fronts on different azimuthal cuts at 171, 193, and 211 A. (f) Same as (e) but for velocities scaled onto 
the lowest cut AO, equivalent to angular velocities in the azimuthal direction with respect to the disk center. 

dip at 171 A and a shoulder at 193 and 211 A, indicating 
a heated sheath surrounding the CME front. 

Within the EUV wave pulse (e.g., panel (n)), it usually 
takes several to 20 minutes to reach the maximum inten- 
sity variation at each wavelength, followed by a com- 
parable or longer recovery lasting up to tens of minutes. 
The maximum 193 A enhancement is about 5-10% above 
the pre-event level, while the maximum 211 A enhance- 
ment and 171 A reduction are about two times larger on 
the 10-20% level and delayed by a few minutes fr om the 
193 A maximum (jLiu et all 120101 : iMa et al.ll20TH ). This 
indicates prompt plasma heating within minutes across 
the characteristic temperatures of the 171, 193, and 
211 A channe ls, from 0.8 MK to 1.6 MK and then 2 MK 
res pectively (IQ'Dwver et al.l 120101 : iBoerner et "all 120111 : 
cf., iKozarev et al.l 1201 lh . The subsequent recovery af- 
ter the maximum toward the pre-event level is a cooling 
phase. Such emission variations were largely reproduced 
in MHD simulations and ascribed to heating by initial 
compression followed by cooling by restoring rarefaction 
in an adiabatic process associated with a compressible 



wave ([Downs et al.ll2012[ ). 

This general trend is evident in the quiet Sun to the 
south of the eruption. To the north, it is complicated by 
organized structures, such as the coronal cavity /filament 
and polar coronal hole boundary producing oscillatory 
intensity variations (panels (d) and (e)). 

B2. KINEMATICS AND GENERAL DECELERATION 

To measure the velocity of the global EUV wave front, 
we visually identified the initial intensity change, regard- 
less of brightening or darkening, at each distance in base- 
ratio space-time diagrams (e.g., Figure 2]) obtained from 
azimuthal cuts shown in FigurelSTb ). We have also tested 
autom ated methods (e.g., see lLiu et al. 2010; L ong et aTl 
l2011b[ ) , but all failed to consistently detect the wave front 
because of the complex intensity variations, as shown in 
Figure [9] Yet, the visual detection, with acceptable un- 
certainties, yielded the most consistent result. 

The resulting EUV wave front positions versus time are 
shown in Figure [TUl for the 193 A channel, as an example, 
and fitted with linear (red dashed) and parabolic (blue 
dashed) functions. In general, both functions fit the data 
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nearly equally well within uncertainties, but parabolic 
fits yield slightly smaller x 2 s and indicate decelerations 
up to to — 0.88 ± 0.17 km s -2 (except for the acceleration 
in panel (a)). At the lowest cut AO (height h = 38"), lin- 
ear fits give an average velocity of 650 ±30 kms -1 for the 
northward wave and ~50% slower at 370 ± 20 kms -1 for 
the southward wave. This north-south asymmetry gener- 
ally holds at different heights. There is a weak increase 
of the velocity with height, but not significant beyond 
uncertainties. 

The EUV wave velocities measured at different wave- 
lengths show substantial scatt er (Figure |Tp[e) and (f)) , 
consistent with those found bv lLong et all (|2008l . l2011af ). 
but their identified trend of higher velocities in cooler 
lines is not present here. At a given height, the velocity 
generally decreases from 193 to 211 and 171 A, likely due 
to the delay of the same trend noted in Appendix Bl. 
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Fig. 11. — 193 A base ratio space-time diagrams from selected 
vertical cuts in Figure Eff), showing height-dependent arrival time 
of the EUV disturbance. Each panel is labeled with the azimuthal 
position s of the vertical cut, as marked by the cross signs in Fig- 
ure [Ifc). The cross signs here on the y-axis mark the heights of 
the azimuthal cuts A0-A3. 
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Fig. 12. — (a)-(n) Running ratio space-time diagrams at 193 A 
from selected radial cuts Ri and vertical cut V6 shown in Fig- 
ures EJa) and (f). Distance is measured from the eruption center 
located at a height of 49" on V6. Note the acceleration of the CME 
front and followup loops at high altitudes and the generation of 
the EUV wave at low altitudes. The red dotted and dashed lines 
are parabolic fits to the CME and EUV wave fronts, respectively, 
labeled with the fitted acceleration (in kms -2 ) and maximum ve- 
locity (in kms -1 , i.e., the initial velocity for deceleration or the 
final velocity at s = 370", top distance of the plot, for accelera- 
tion). The vertical white dotted line marks the onset of the rapid 
CME expansion and the flare impulsive phase at 23:05 UT. (o) The 
bottom panel shows the angular distribution of the CME acceler- 
ation measured from these cuts. The angle is measured from the 
local vertical in the counter-clockwise direction. 



B3. HEIGHT DEPENDENCE AND FORWARDLY INCLINED 
WAVE FRONTS 



Quasi-periodic Wave Trains Within A Global EUV Wave Traveling Through A Coronal Cavity 



15 



The height dependence of the EUV wave can be exam- 
ined through its arrival at vertical cuts defined in Fig- 
ure |3](f). As shown in the resulting space-time diagrams 
of Figure [HJ not far from the eruption (top three rows), 
a general pattern, similar across different AIA channels, 
is as follows. In the low corona (h < 300" = 220 Mm), 
the initial intensity variation forms a stripe of a nega- 
tive slope, indicating that the arrival of the disturbance 
at the vertical cut first occurs at a high altitude and 
then progresses toward lower altitudes; in the high corona 
(h > 300"), the trend is opposite. For example, over a 
height range of 200", the downward progression at cut V3 
(south of the eruption) results in a 3 min delay, trans- 
lating to an apparent vertical velocity of — 760 kms -1 , 
while the upward progression at cut V7 (north of the 
eruption) produces a 4.5 min delay and an apparent ve- 
locity of 540 kins" 1 . This pattern means that the leading 
front of th e EUV wave near th e eruption has a sphere- 
like shape ([Veronig et al.ll2010f ): viewed from the erup- 
tion center, the front is upwardly concave in the high 
corona and forwardly inclined toward the solar surface 
in the low corona, as further discussed in Section [6] 

The delayed arrival at the vertical cut V8 shown in 
Figure [5fk), for example, can be used to reconstruct the 
corresponding wave front. For each height h, we first in- 
terpolated the azimuthal wave velocity v(h) from those 
average velocities measured on cuts A0-A3 at different 
heights (see Figure [TO]). We then obtained the azimuthal 
distance traveled from the actual arrival time t(h) to a 
reference arrival time of t(0) =23:16:55 UT at the coro- 
nal base (h = 0), As(h) = [t(0) - t(h)]v(h). This dis- 
tance was added to the azimuthal position s cut (h) of this 
vertical cut measured from the eruption center, to ob- 
tain the position of the instantaneous wave front at £(0), 
s(h) = s cut (h) + As (ft). Repeating this for all heights 
and then transforming the azimuthal positions back to 
the Cartesian coordinates gave the resulting wave front, 
as shown in red in Figure 03(1). 

APPENDIX C 
CME KINEMATICS AND EARLY EXPANSION 

Figure [12] shows the expansion in different directions 
of coronal loops involved in the CME, as captured in 
radial cuts defined in Figure [3ja). Strong acceleration 
(a > 0.4 km s -2 ) only occurs within an angular range 
—50° < < 70° of the local vertical at the eruption 
center (panel (o)). In this range, the CME experiences 
a two -phase, slow to fast expansion ([Sterling fc Moorel 
|2005[ ). the transition of which occurs at 23:05 UT, the 
beginning of the flare impulsive phase. The early-phase 
has an acceleration <0.1kms -2 and maximum veloc- 
ity ~100kms -1 (panel (a)). The largest acceleration of 
0.70±0.06 kms -2 in the late-phase takes place in the ver- 
tical direction (panel (h)), reaching 240±20kms _1 at the 
edge of the AIA FOV. Comparable acceleration possibly 
continues into the SOHO/LASCO FOV (1.5-30 #©), as 
shown by the dotted line in the inset of Figure EJc). The 
CME then decelerates at a = -0.019 ±0.005 km s~ 2 with 
an initial and average velocities of v max = 990±70 kms" 1 
and (v) = 820 ±40 kms -1 , according to the CME catalog 
(https: / / cdaw.gsfc.nasa.gov/CMEJist) . 

Beyond the angular range of -50° < < 70°, the CME 
acceleration remains below 0.3 kms -2 , with no two-phase 
distinction but a north-south asymmetry. To the north 



of the vertical, the expanding loops gradually slow down 
and/or fade away (panel (e)). This first occurs at the 
lowest cut and progresses toward cuts at greater heights, 
as indicated by the short yellow arrows. To the south, the 
early slow expansion is absent, but a sudden brightening 
caused by a rapid downward expansion of some loops 
(see Movie 1C), indicated by the short orange arrows, 
takes place sequentially from cut R4 to R0. It is at the 
outermost edge of this brightening near the coronal base, 
so = 150" = 110 Mm from the eruption center, that the 
southward global EUV wave is launched at 23:14:20 UT 
(panels (m) and (n)). Its northward counterpart is first 
detected at a similar distance 3 min earlier. 

We note in panel (h) that, in the vertical direction, 
some loops at the eruption center are pushed downward 
from an initial height of fto = 150" (with the 49" height 
of the origin of the cut included), followed by transverse 
oscillations. In lateral directions, at the onset of the rapid 
CME expansion (vertical dotted line), all participating 
loops are also located approximately beyond a distance 
so = 110 Mm from the eruption center. This indicates 
that the rapid expansion involves the entire active region 
of tar = 110 Mm in radius and originates at its center 
of height fto = 110 Mm. These observations imply that 
this expansion, with a significant downward component, 
may play a key role in driving the EUV wave (see more 
discussion in Section |CT2]) . 

APPENDIX D 
ESTIMATES OF EUV WAVE ENERGY AND CAVITY 
MAGNETIC FIELD 

We can utilize the oscillations of the coronal cavity 
to estim ate the energy bud get of the global EUV wave 
(see, e.g., Bal lai et al .120051 ) and the cavity magnetic field 
strength, i.e., by coronal seismology. These are order-of- 
magnitude estimates with comparable uncertainties. 

The cavity has a diameter D cav = 135" = 99 Mm mea- 
sured on azimuthal cut A2 in Figure [5](c). Its embed- 
ded filament, likely underneath the axis of the flux rope 
runn ing along its full length (e.g., lLow fc Hundhausenl 
Il995[ ). has been quite stable for a week before this 
event. It was best seen a few days ear lier (say, Septem- 
ber 4, at http: / / www, solar monitor .or g[ ) , located at N50° 
on the disk, spanning a longitude range of 35° and a 
length of L cav = 273 Mm. Assuming the cavity being 
a cylinder of such a dimension, its volume is V cav = 
7r(D cav /2) 2 L cav . We adopt an average electron density 
of n e = 1 x 10 9 cm -3 of the tenuous cavit y and its 
dense filament combined (jVasquez et al. I [20091 ) and thus 
p — n e m p — 1.7 x 10 -15 gcm -3 . Its total mass is m cav = 
pVcnv = 3.5xl0 15 g com parable to a typical CME mass 
(jAschwanden et alTl2QQ9h . Now, take a velocity ampli- 
tude of v osc = 9 kms -1 , we obtain the kinetic energy of 
the cavity of £^, C av = m C avvl sc /2 ~ 2 x 10 27 ergs and en- 
ergy per unit length along its axis of e kjCav = E kjCav /L cav . 
We further assume such an energy is uniformly dis- 
tributed in space at the same projected distance from 
the eruption center, s = 450" = 329 Mm, measured on 
cut A2 (Figure 03(c)), as a consequence of a surface-like 
wave in the low corona. Then, the global EUV wave en- 
ergy is E wv = efc >cav (27rs) ~lx 10 28 ergs. Both energies 
are 3-4 orders of magnitude smaller than the kinetic en- 
ergy Ek y cME ~ 1 x 10 31 ergs of the CME, assumed at the 
same mass as the cavity moving at 820 kms -1 measured 
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by SOHO/LASCO. This wave energy is, however, 100 
times greater tha n that of another EUV wave and of a 
typical nanoflare (B allai et al.ll2005f ). 

These transverse oscillations are likely fast kink modes. 
In this case, the period for the Nth harmonics is P = 
2Z/(7Vc/e), where I is the loop length and c& ~ va ~ Vf is 
the kink speed in \ow-f3 plasma with a loop density of the 
same order of the surrounding plasma (va and Vf being 
the Alfven and fast-mode speeds, respectively). Using 



the measured average period of 27.6 min for the fun- 
damental mode N = 1 and assuming c& = 980 km s -1 , 
we have I = Pck/2 w 800 Mm w 3L cav , where L cav = 
273 Mm is the cavity's axial length measured above. This 
gives a pitch angle of = arccos(L cav //) = arccos(l/3) ~ 
70°, implying an expe cted, highly twisted flux rope 
([Carter fc Erd elyi 2008). If we adopt a typical density of 



n e = 2 x 10 8 cm (p = n e m p = 3 x 10 g cm" , from 
iSchmit fc Gibsonll2011[ ) for the coronal cavity, we obtain 
a magnetic field strength of B = va y/knp ~ 6 G. 
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